Self-diffusion of Rod-like Viruses Through Smectic Layers 
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We report the direct visualization at the scale of single particles of mass transport between 
smectic layers, also called permeation, in a suspension of rod-like viruses. Self-diffusion takes place 
preferentially in the direction normal to the smectic layers, and occurs by quasi- quantized steps of 
one rod length. The diffusion rate corresponds with the rate calculated from the diffusion in the 
nematic state with a lamellar periodic ordering potential that is obtained experimentally. 



Since the pioneering work of Onsager on the entropy 
driven phase transition to a liquid crystalline state 
the structure and the phase behavior of complex fluids 
containing anisotropic particles with hard core interac- 
tions has been a subject of considerable interest, both 
theoretically [2] and experimentally [3]. Understanding 
of the particle mobility in the different liquid crystalline 
phases is more recent [4] . In experiments various methods 
have been applied to obtain the ensemble averaged self- 
diffusion coefficients in thermotropic [5[ and amphiphilic 
[q] liquid crystals, block copolymer pT^and colloidal sys- 
tems 0. Only a few studies have been done where dy- 
namical phenomena are probed at the scale of a single 
anisotropic particle: the Brownian motion of an isolated 
colloidal ellipsoid in confined geometry [9] and the self- 
diffusion in a nematic phase formed by rod-like viruses 
[lO] represent two recent examples. In the latter case, the 
diffusion parallel {D\\) and perpendicular {D^) to the av- 
erage rod orientation (the director) has been measured, 
showing an increase of the ratio D\\/D^ with particle con- 
centration. Knowledge of the dynamics at the single par- 
ticle level is fundamental for understanding the physics of 
mesophases with spatial order like the smectic (lamellar) 
phase of rod-like particles. In this mesophase the particle 
density is periodic in one dimension parallel to the long 
axis of the rods, while the interparticle correlations per- 
pendicular to this axis are short-ranged (fluid-like order). 
For parallel diffusion to take place, the rods need to jump 
between adjacent smectic layers, overcoming an energy 
barrier related to the smectic order parameter [11]. This 
process of interlayer diffusion, or permeation^ was first 
predicted by Helfrich [12]. In this Letter, we use video 
fluorescence microscopy to monitor the dynamics of indi- 
vidual labeled colloidal rods in the background of a smec- 
tic mesophase formed by identical but unlabeled rods. In 
this way we directly observe permeation of single rods in 
adjacent layers. As in the nematic phase, self-diffusion in 
a smectic phase is anisotropic: the diffusion through the 
smectic layers is shown here to be much faster than the 
diffusion within each liquid-like layer, i.e. D\\/D^ ^ 1, 
in contrast to thermotropic systems. Moreover, since the 
individual rod positions within the layer are monitored, 
the potential barrier for permeation is straightly deter- 
mined for the first time. The permeation can then be 
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FIG. 1: (a) Time sequence of an overlay of fluorescence and 
DIG images showing labeled particles jumping between ad- 
jacent smectic layers (At = 0.071s is the time between two 
frames). The layer spacing is L 0.9 /xm. (b) Displacement 
of a given particle in the direction parallel (red) and perpen- 
dicular (black) to the director. The green lines indicate the 
residence time, i.e. the time for which one particle stays in a 
given layer. 



described in terms of Brownian particles diffusing in a 
one-dimensional periodic symmetric potential. 

The system of rods used in this work consists of fila- 
mentous bacteriophages fd^ which are semi-rigid polyelec- 
trolytes with a contour length of 0.88 jam, a diameter of 
6.6 nm, and a persistence length of 2.2 /j^m 0]. Suspen- 
sions of fd rods in aqueous solution form several lyotropic 
liquid crystalline phases, in particular the chiral nematic 
(cholesteric) phase and the smectic phase 



IJ. The ex- 



istence of a smectic phase in suspensions of hard rods is 
an evidence of the high monodispersity and therefore of 
the model system character of such filamentous viruses 
[l^ . The colloidal scale of the fd bacteriophage fa- 
cilitates the imaging of individual rods by fluorescence 
microscopy, as well as smectic layers by differential inter- 
ference contrast (DIG) microscopy [l|. Fig. IHa) shows 
a sequence of images of a single region [17] where both 
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FIG. 2: Probability density function in space at different 
times for two ionic strengths: 20 mM (a) and 110 mM (b). 
The functions are normahzed to one, the z-axis is scaled by 
the smectic layer thickness L. 



techniques are combined. A comparison of the images 
shows that some rods jump between two layers while oth- 
ers remain within a given layer. The trajectory of one of 
the rods is plotted in Fig. [TJb) in the direction paral- 
lel (z) and perpendicular (x) to the director. This figure 
summarizes the key observation of this Letter: the diffu- 
sion throughout the smectic layers takes place in quasi- 
quantized steps of one rod length i.e. the mass transport 
between the layers is a discontinuous process. Moreover, 
it shows that the diffusion within the smectic player is 
extremely slow [18j. 

The "hopping-type" diffusion is mainly the conse- 
quence of the underlying ordering potential of the smec- 
tic phase and the vacancies available in adjacent layers. 
A phenomenological expression for permeation has been 
derived by coupling the displacement of a segment of a 
smectic layer u to the compressibility modulus B via the 
permeation parameter [llj : 



du 
'dt 



(1) 



On a single-particle level, the fundamental solution of 
this diffusion equation is the self- van Hove function [3], 
which is the probability for a displacement z during a 
time t: 
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Since single particles are experimentally identified, the 
self-van Hove function can be directly obtained from the 
histogram of particle positions after a time t, as plotted in 
Fig. Elfor low (I = 20 mM) and high ionic (I = 110 mM) 
strengths. For a fluid made of Brownian particles, a 
smooth gaussian distribution that smears out over time 
is expected for the self-van Hove function. However at 
low ionic strength, G{z^ t) shows distinct peaks exactly at 
integer multiples of the particle length (and therefore of 
the layer thickness, see Fig. [2fa)) as expected from visual 
observation (Fig. [1]). At high ionic strength the curves 
are smoother (Fig. EJb)), but in all cases the experimen- 
tal self- van Hove function is not gaussian at any time. 
This implies that the permeation parameter A5 in Eq. [T] 
is a function of position z^ due to the energy landscape 
imposed by the smectic layers. 

The energy landscape can be determined experimen- 
tally from the distribution of particle positions with re- 
spect to the middle of a layer parallel to the director. To 
this end, time windows are selected where the particle 
remains for ten frames or more within the same layer. 
The distribution of particles within a single layer is then 
obtained by addition of all particle positions relative to 
the average position of particles for all selected time win- 
dows. The resulting distributions are plotted in Fig. [Sja) 
for the two ionic strengths. To obtain the total particle 
distribution for the full smectic phase, the distributions 
of particles in a single layer (Fig. [3ja)) is added period- 
ically to itself at all integer numbers of layer spacing L 
(Fig. [ifa)). The smectic ordering potential is then de- 
duced from the Boltzmann factor P{z) ^ ^-Uiayer{z)/kBT 
for the probability of finding a particle at position as 
shown in Fig. [Sjb). Both potentials can be best fitted 
with a sinusoidal Uiayer{z) = /7o sin(27rz/L), giving an 
amplitude of Uq = 1.36 ksT at low ionic strength and 
Uq = 0.66 ksT at high ionic strength. The difference 
between the two amplitudes explains the fact that for 
I = 20 mM the self-van Hove function exhibits discrete 
peaks, while for I = 110 mM the potential barrier is small 
enough to exhibit a monotonic behavior of the probabil- 
ity density function. The reason for the more pronounced 
potential at low ionic strength might be that electrostatic 
interactions between rods are more long-ranged, i.e., par- 
ticles are more strongly correlated so that it is more dif- 
ficult to create a void between them. The fact that the 
potential can be fitted by a sinusoidal is remarkable by 
itself. Indeed, the use of such a potential is very common 
due to its simplicity [20] , but this ordering potential has 
never been directly observed until now. Moreover the 
height of the potential, i.e. the smectic order parameter, 
can be directly obtained. 
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FIG. 3: (a) Histogram of time averaged particle positions par- 
allel to the director within the smectic layer at 20 mM (■) 
and 110 mM (o). The histogram is normalized by the total 
number of positions, (b) Resulting effective mean ordering 
potential in the z-direction obtained by applying the Boltz- 
mann factor. The solid lines are a fit to a sinusoidal potential. 



The overall mean square displacement (MSD) of rods 
parallel and perpendicular to the director of the smec- 
tic and nematic phase is plotted in Fig. [4] for both 
ionic strengths. The time evolution of the MSD given by 
(Ar^(t)) rsj provides the diffusion exponent 7: 7 < 1 
is characteristic of a 5iiWiffusive behavior, while 7 > 1 is 
referred to as 52/pe?xliffusion. The parallel motion is close 
to be diffusive (7 1) in both the nematic (7 = 0.97) 
and smectic (7 = 0.94) phases for I = 110 mM and in the 
nematic phase for I = 20 mM (7 = 0.95). Only the paral- 
lel motion in the smectic phase for low ionic strength, i.e. 
where the discrete peaks in the self- van Hove function are 
observed, is significantly subdiffusive: 7 = 0.81. The per- 
pendicular motion is in all cases strongly subdiffusive: for 
I = 110 mM, 7 reduces from 0.63 before to 0.56 after the 
nematic-smectic (N-Sm) transition and for I = 20 mM it 
reduces from 0.68 to 0.46. Anomalous subdiffusive be- 
havior has often been observed in systems where diffu- 
sion takes place by steps, e.g. in case of release from a 
surrounding cage [21]. This "cage escape" might be at 
the origin of the observed subdiffusive behavior for both 
parallel and perpendicular diffusion. For parallel diffu- 
sion the cage is formed by the energy barrier imposed by 
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FIG. 4: (a) Mean square displacement parallel to the director 
for the two indicated ionic strengths in the nematic phase 
(red line), in the smectic phase (blue dashed line), and in 
the nematic phase considering the oscillatory potential (red 
dotted line), (b) Mean square displacement perpendicular to 
the director (same convention as above) . The insets show the 
same data in a log-log scale, yielding the degree of subdiffusion 
from the linear regression (black line). 



the smectic layers, as shown by smaller 7 for higher or- 
dering potential. Perpendicular diffusion at high volume 
fractions is only possible through a reptation-like motion 
along the long axis to escape the local excluded volume, 
as observed for polymers for which typically 7 = 0.5 [22[. 
This excluded volume is huge, even for thin rods at high 
orientational order, due to the large rod aspect ratio of 

130. In addition, perpendicular diffusion in the smec- 
tic phase is hindered due to the ordering potential, which 
couples this diffusion to the permeation and which thus 
explains the decrease of 7 from the nematic to the smectic 
phases. For subdiffusive systems, a non- Gaussian distri- 
bution of the probability density functions has been ob- 
served as in Fig. [2] [2T], even though these two features 
are not a priori correlated. Note also that boundary ef- 
fects might influence the probability density [23]. 

The anisotropy in the total diffusion, D\\/D_\_, which 
is about 20 in the nematic phase [10], increases in the 
smectic phase as a result of the pronounced subdiffusiv- 
ity of the perpendicular motion (decrease of 7). These 
observations show an opposite trend as compared to ther- 
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inotropic liquid crystals 0, [1], where usually D\\/D^ 
evolves from being larger than one at temperatures close 
to the N-Sm transition temperature to being smaller than 
one at lower temperatures [24] . Therefore the diffusion in 
the smectic phase can be effectively considered as a one- 
dimensional diffusion of a Brownian particle in a periodic 
potential in the high friction limit. A general expression 
for such a diffusion process is given by [25,]: 



The brackets indicate averaging over one period of the 
ordering potential. The diffusion coefficient in the smec- 
tic phase can then be calculated taking Dq as the diffu- 
sion coefficient in the nematic phase close to the N-Sm 
transition, and using Uiayer as obtained from the fit of 
the potentials plotted in Fig. [3l the diffusion coefficient 
decreases by a factor 0.84 at I = 110 mM and by a fac- 
tor 0.44 at I = 20 mM. Indeed the MSD in the smectic 
phase is obtained from the MSD in the nematic phase, 
using these factors for both ionic strengths (see Fig. [4]), 
although at I = 20 mM some deviation appears due to the 
subdiffusivity in the MSD. Thus, we have shown how the 
mobility of rods decreases after the N-Sm transition, con- 
trary to the isotropic-nematic transition where the global 
mobility increases due to entropic gain [l|, [l^ . It seems 
therefore to indicate that fd virus suspensions do not be- 
have as a system of rigid hard rods for high concentration 
in agreement with a recent work [l^. Moreover, the very 
slow diffusion within the layers suggests that the smectic 
phase of semi-flexible colloidal rods consists of layers of 
glass-like, rather than fluid-like, particles. 

In conclusion, we have for the first time visualized the 
process of permeation in the smectic phase at the scale 
of single particles for a system of charged rods. This 
allowed us to give a full and coherent description of the 
diffusion process without any assumptions on the system. 
The diffusion is strongly anisotropic in the direction nor- 
mal to the smectic layers and quasi-discontinuous due to 
the presence of the layers. The parallel diffusion rate 
complies with the rate in the nematic phase, taking into 
account the ordering potential, which is obtained directly 
from our measurements. 
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